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Abstract 

Consider n + m non- intersecting Brownian bridges, with n of them 
leaving from at time t = —1 and returning to at time t = 1, while 
the m remaining ones (wanderers) go from m points to m points 
5j. First we keep m fixed and we scale ai,bi appropriately with n. 
In the large-n limit we obtain a new Airy process with wanderers, in 
the neighborhood of -v/2n, the approximate location of the rightmost 
particle in the absence of wanderers. This new process is governed by 
an Airy-type kernel, with a rational perturbation. 

Letting the number m of wanderers tend to infinity as well, leads 
to two Pearcey processes about two cusps, a closing and an opening 
cusp, the location of the tips being related by an elliptic curve. Upon 
tuning the starting and target points, one can let the two tips of the 
cusps grow very close; this leads to a new process, which we conjecture 
to be governed by a kernel, represented as a double integral involving 
the exponential of a quintic polynomial in the integration variables. 
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1 Introduction 

Consider n + m non-intersecting Brownian motions (Brownian bridges) on 
R depending on time t G [—1, 1], with n of them leaving from and returning 
to 0, while the m remaining ones leave from < • • • < oi and are forced 
to end up at 6m < • • • < We denote by Xi{t) the position at time t of the 
ith largest Brownian particle among the n + m non-intersecting Brownian 
bridges. Denote by T) the conditioning event defined by the following condi- 
tions: 

(i) non-intersecting paths: Xi{t) > X2{t) > . . . > Xm+n{t), t & (—1, 1), 

(ii) n bridges from to 0: 1) = Xj(l) = for i = m + 1, . . . ,m + n, 

(iii) m wanderers from to bf. Xi{—1) = ai, Xi{l) = bi for i = 1, ... ,m. 
Then, denote the conditional probability under V by Pab, i-e., 

F,^{-)=F{-\V). (1.1) 

The interest in non-intersecting Brownian motions stems from a paper 
by Dyson [19], who made the important observation that putting dynam- 
ics into the GUE-random matrix model (Ornstein-Uhlenbeck Processes on 
the real and imaginary parts of the entries) leads to finitely many non- 
intersecting Brownian motions on M for the eigenvalues (stationary process). 
A space-time transformation enables one to map the above Dyson process 
into non- intersecting Brownian motions starting from and returning to 0; 
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Figure 1: Non- intersecting Brownian bridges with m wanderers, leaving from 
a = Xq(1 + dn~^^^) and forced to 6 = Xq(1 — bn~^^^), with a < 6, where 

X- = v^y|±fa, x^ = V2n^J^^. The dotted hne hnking {xq,-1) to 
{xq, 1) is tangent to the curve x = ^/2n{l — t"^) at the point {xq, to)- 



see formula (1.7) in [1]. In their work on coincidence probabilities, Karlin- 
McGregor [30] found a determinantal formula for the transition probabil- 
ity of non-intersecting Brownian motions. The relationship between non- 
intersecting Brownian motions, matrix models and random matrix theory 
has been developped starting with Johansson [26] and has led to the Airy 
and other processes [2-4, 17,38,39], when the number of particles tend to 
infinity, see also [18]. 

At first, consider the motion of the non- intersecting Brownian particles 
above, but with m = 0, and let n become very large. The Airy process A{t) 
describes this cloud of particles ("infinite-dimensional diffusion"), but viewed 
from any point on the "edge" C : x — >/2n(l — t^) of the set of particles, 
with time and space properly rescaled; the Airy process will be independent 
of the point chosen and will be governed by the Airy kernel. This process 
was found by Prahofer and Spohn [36] in the context of stochastic growth 
models and further investigated in [2,22,27,28,38]. 

Assume now a fixed and finite m > 1 and all ai — 0, with the target 
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points all equal to b scaled as 6 = po\'2n > 0. Does it affect the Brownian 
fluctuations along the curve C for large n? No new process appears as long as 
one considers points {y, t) G C, below the point of tangency of the tangent to 
the curve passing through (po\/2n, 1). At this tangency point the fluctuations 
obey a new statistics, which we call the Airy process with m outliers ^^(r), 
governed by a rational perturbation of the Airy kernel, see [1]. This kernel 
was already considered by Baik-Ben Arous-Peche [5] and Peche [35] in the 
context of multivariate statistics. 

The first result in this paper concerns the limiting process, described in 
(11.11) . in the large-n limit, while keeping m fixed; this process is denoted by 

At''\r). This paper deals with the statistical fluctuations of the edge of 
the cloud of particles near any point on the curve m 
the presence of wanderers. To do so, consider the tangent line to the curve 
C, with point of tangency (0:0,^0)5 as in Figure [U this tangent intersects 

the lines t = —1 and t = 1 at the points Xq = = \/2nJ^^ and 



~ ~ v^2ny Y+l^ respectively. Consider now m wanderers leaving 
from neighboring points (when n gets large) of the point Xq at time t = —1 
and forced to neighboring points of Xq at time t = 1. A first part of this 
paper is to show that the fluctuations near the edge of the cloud and near the 
point (xo,to) obeys a new statistic, independent of the point (xo,to) chosen 
on the curve above, showing universality within that class. 

At a first stage (Theorems 11.11 and II. 2p . the result will be shown for a 
vertical line tangent to C at the point (-\/2n, 0), whereas Theorem 11.31 deals 
with the universality result. The non-intersecting nature of the first n bridges 
implies that the largest one will again reach a height of about \/2n. So, it is 
natural to consider the following scaling of the starting and the target points 

2n(l + ^J and 6. = V2n |^1 - | . (1.2) 

With this scaling, the m wanderers will interact with the bulk (of n particles, 
with n very large) in a non-trivial way, upon considering regions close to 
X = \/2n and t = 0, namely at space-time positions (x, t) which scale like 

t = rn-^l^^ X = V2^ + , . (1.3) 

This will only be so under some geometric condition: the lines connecting 
the starting and target points in {x, t)-space must pass to the left of \/2n 
at t = 0; see Figure [H Then, the first result concerns the gap probability 
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at a given time r for very large n and keeping m finite and fixed, i.e., the 
probability that a set is not visited by any of the n + m Brownian bridges 
at time r. Thus, in Theorems 11.11 and 11.21 a different (non-trivial) process 

A^^''\t) will appear due to the interaction of the m wanderers with the Airy 
field in the neighborhood of (x,t) = (\/2n, 0). Note that in the absence 
of wanderers the particles must look, near the edge, like the Airy process. 
This also explains why the kernel fll.Sp obtained below is another rational 
extension of the Airy kernel. 

Theorem 1.1. Consider points ai and hi, as in with (im < ■ ■ ■ < di < 

bi < . . . < bm on the real lin^. Given any compact set C M, the gap 
probability at rescaled time-space U.3\) is given, in the large-n limit, by 

Jim X, (^) G + ^^}) = det(l - XeK'Jxe)^^ 

= P(^(^'^)(r)nE = 0), (1.4) 

where Xe{0 = ^ E), where det denotes the Fredholm determinant on 
L^(R) and where the kernel K'^ is given by 



(1.5) 

The integration contours are as follows: Ta> goes from e ^'^'/^oo to e^'^'/^oo, 
and passes on the right of all the ai — t, while T^^ goes from e'^'^^^oo to 
g-TTj/s^Q^ anc? passes to the left of all b^ — t. Moreover, the two contours do 
not intersect; see Figure\^for an illustration. 

This kernel has also appeared in recent work of Borodin and Peche [15], 
as a limit of a directed percolation in a quadrant with defective rows and 
columns, itself a generalization of a kernel of Baik-Ben Arous-Peche [5] and 
Peche [35] and considered in [1] in the context of non-intersecting Brownian 
motions. The same limit process occurs in the asymmetric exclusion process, 
see [14,23]. The proof of Theorem 11.11 will be given in Section [21 when the 
points Oj and the points bi are all different. When the a^'s all coincide, and 
similarly the 6j's, the proof of Theorem 11.11 breaks down and must be replaced 
by another one; two approaches are being discussed here (see Section H]): (1) 
using a certain moment matrix, (2) using biorthogonal polynomials. 



^The inequalities that all the ai be smaller than all the hiS means geometrically that 
the lines connecting corresponding points intersect the a;-axis to the left of a: = V2n; see 
Figure [D 
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Figure 2: Integration paths of the kernel of (11 .Sp . 



In Theorem 11.21 (see Section [3]) the first result will be extended to the 
joint gap probabilities at different (rescaled) times ri, . . . , r^. Obviously The- 
orem 11.11 is the specialization of Theorem 11.21 to the one-time case. 

Theorem 1.2. Consider i distinct times ti,T2, . . . ,Ti and compact sets 
Ei,...,EiCR. Then, 



\A;=1 



det(l - XEK':i>XE) = P ( n 'H^'^) n = 0} ) , (1.6) 



\k=l 



where XEijkiC) '■= ^ Ek)- Here det denotes the (matrix) Fredholm 
determinant on the space L^({ri, . . . r^} x M) and the extended kernel is 
given by 



A/47r(r2 - Ti) 

^1 ^ 'z:rf'i.^>}^^ ■ 



The integration contours are as in FigurelE, but with ak — r replaced by a^ — r2 
and hk — t replaced by 6^ — ri . 

A similar statement can then be made along any point (xo,to) of the 
curve X = A/2n(l — t^), with tangent intersecting the lines t = — 1 and t = 1 
at the points 



Xn = = V 2n\ and Xn = = v 2n\ , 1.8) 



respectively. 

Theorem 1.3. (Universality statement) As before, consider i distinct 
times Ti, T2, . . . ,T£ and compact sets Ei, . . . , Eg <Z M.. Also, consider m Brow- 
nian wanderers, now leaving from the points ae = Xq (1 + and forced to 

bi = Xq (l — ^Jts), with the conditio"^ 0,^ ^ ■ ■ ■ ^ 0,1 < bi < . . . < bm- For 
n large, pick i points in a rT^I^ -neighborhood of {xQ.to), lying on the curve 

Xk := ^2n{l-tl), with tk := to + -^^^ItI^, 1 < k < i. (1.9) 
Then, the following limit hold^: 




(1.10) 



Remark 1.4. For (xq, to) = ("\/2n, 0), this statement reduces to Theorem ll.2l 
as can seen from Footnoted 

In view of the new process Am''\T), it seems natural to let the number 
of wanderers m to go infinity. For simplicity, consider the case where the m 
wanderers all start from the same point a, and end up at the same point b, 
with a < b, with the scaling 

a = am^/^, b = /5m^/^ with a < p. (1.11) 

Under this scaling, the set of m wanderers itself produces an Airy field which 
then interact with the one already present after the n 00 limit. Thus, 
we might expect that there will be two regions where the Pearcey process 
arises. Indeed, the first Pearcey process occurs when the "Pearcey cusp" 
closes, while the second does when the cusp opens, as illustrated in Figure El 

^Here also, the inequalities that all the di be smaller than all the 6i's means geometri- 
cally that the lines connecting corresponding points intersect the horizontal line through 
{xo,t()) to the left of {xo,to); see Figure [Tj 

■^Expanded out, (^1 + 2^^) reads 
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n — *• cxo 



m ^ oo 



Figure 3: Illustration of the two Pearcey processes, arising around the two 
cusps. 




Figure 4: Integration paths of the Pearcey kernel K]^ defined in f ll.l2p . The 
two solid lines form together Fx, the dashed line is the ^-integration path. 



The reader is reminded of the extended Pearcey kernel K^{6i, vi, 62, V2) 
with space-time parameters {6i,Vi), which is given by 



1(^2 > 



+ 



271(62 - Oi] 

pioo 



{2m) 



dz / dz- 



I g-54/4+eiz2/2-j;ii ' 



'1.121 



where the path Fx is illustrated in Figure HI see Tracy- Widom [39]. This 
leads to Theorem 11.51 established in Section 
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Theorem 1.5. Let the starting point a and the target point b of the m wan- 
derers for the Airy process with wanderers /[1.1(J\) grow with m, as a = am}/^ 
and h = f3m^^^ with arbitrary a < (3. Given a < (3, the following equations, 

(3 — a = , with (x,a) G £ : Aa^x^ — 2x + 3 = 0, (elliptic curve) (1.13) 

2 — X 

have a 

unique solution (x, a) := (x,cr+) G ((§,2) x (— 1,0)) (opening cusp) 
unique solution (x, cr) := (x, (T_) G ((§7 2) x (0, |)) (closing cusp). 

Then, the Airy process with m wanderers J^m^\T) properly rescaled as 
m oo, converges to two (identical) Pearcey processes V {9) about two cusps, 
one opening cusp (T^) and one closing cusp (T^) about^ 

Tr~.T^m^l\ ir^Xm^'\ ujzth T^:=^^-^, X:=4(l-2x), 
± , s , ^2 2-x ^ 

(1.14) 

with < ^ < T+. To be precise, upon using the two different scalings 
Iil.l6\) below, depending on the opening or closing cusp, one has, for any 
i = 1,2,..., that the limit of the gap probability of the sets Ei, . . . ,E£ at 
times Ti, . . . ,Ti is given by the same (matrix) Fredholm determinant, 



lim F{n\A^r''''^"''''\n)nEk 

n. — >r\n \ ' 'I 



det (1 - XeK^Xe),..,^_,.,^. =: P n i'^i^^) H = 0} , (1.15) 



where the rescaling from the space-time variables {Ei,Ti) to the new space- 
time variables {E^, 9i) is imposed by the initial scaling to yield 



1 o„ /2(x-l)^'/' 



2 V l^±|a;' / (1-16) 

Ei = Xm2/3 - /tVi^im^/^ - nEim-^/^^ 

Remark 1.6. Note that the involution: vi V2, 9i ^ —92, <^ T_, 
(T+ ^ (T_ = —0-+, where Vk G E^, maps the opening cusp into the closing 
cusp and, in particular, acts on the kernel fll.l2p to produce the kernel going 
with the closing cusp. 
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n — >■ oo m — oo 

Figure 5: When two Pearcey cusps touch, there will be a new process. 



The tips of the two cusps in Theorem 11.51 come together, when a, (3 —>■ 0, 
and hence x —>■ 3/2, cr± ^ and T± —>■ 0; this is not the only way for this to 
happen, as will be mentioned below. At the very point where the two cusps 
meet, a new process will emerge (as in Figure [5l), which we conjecture to be 
governed by a "quintic kernel". 

Conjecture 1.7. The gap probability for the new process appearing in Fig- 
ure\E is given by the Fredholm determinant of the following quintic kernel.- 

K^(9,r];x,y) = , dz dz r „-5,^ n-3,n — ^jt^, (1-17) 

where the z and z-integration paths are given by the z and z-paths in Figurel^ 
with the orientation indicated. 

To give some evidence to this conjecture, we first notice that the curve S 
(introduced in (11.131) ) contains another real point (besides the real segments 
introduced just after (I1.13P ) namely at (x, cr) = (oo,0), for which (a,/3) = 
(2^/^,-2^/^); there the critical point Wc of the associated steepest-descent 
F-function becomes order 5, with (X, T) = (—2^/^, 0), rather than order 4 as 
in the Pearcey case; this expresses the fact that the two tips come together. 
For this choice of (a,/3) = (2^/^ -2^/^), the source and the target points 

a = V2n ^1 H — ' ^ ~ V2n (l ^ j , with a = am^^"^ and b = [3m} 



corresponds to (7+ < and T_ corresponds to cr_ > 0, with obviously 
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/ 

/ 

/ 

/ 

/ 

4 

I 

I 

Figure 6: Integration path C (dashed hne) and C (sohd hne) of the quintic 
kerneL 



do not, of course, satisfy the inequahty d < b, but rather the opposite in- 
equahty. We then perform an analytic continuation of the (one-time) kerneH 



1 /■ r c—^ /3+£,2i^ — - ~ — f 

{2mf Jt-^^ Jt^^ V3+Ci- - uj 

by moving d and h in the complex plane from their original position d < h 
to a new position 6 < a on the real line. Then by picking d = am}/^ and 
h = l3m}/^^ with (a,/?) = (2-*^/^, —2^/'^) and letting m oo, we show the 
kernel f ll.lSp tends to the quintic kernel (11.171) with the precise contour of 
integration in the figure above. Some evidence in favor of this conjecture 
is given in Section [HI which contains two rigorous statements, with proofs. 
However this does not suffice to prove the conjecture; e.g., it is still unknown 
whether the Fredholm determinant of the quintic kernel fll.l7p determines a 
probability. For numerical methods, see, for instance, Bornemann [9]. 

Remark 1.8. It is interesting to put the three kernels in parallel. Airy, 
Pearcey and quintic together with their appropriate contours, as in Figure [71 



'In the one-time case, one can just absorb the time r in the S, and b. 
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7r/3 /TTx 7r/3 
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Figure 7: Integration paths C (dashed hne) and C (sohd hne) of the (a) Airy 
kernel X-^, (b) Pearcey kernel ^ and (c) Quintic kernel K^. 



i^'^(Ti,Ci;T2,6) = t:-— 2 dz dz 



(27rz)2 Jc z- 



{2mf Jc Jc + T2) -{z + Ti) e^'/3-«2-' 

y^47r(T2 - Ti) 



y^2n{T2 -Ti) 
Iff 1 g225/5_rz3/3_^22^^^_j. 
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2 Airy process with wanderers leaving from 
and going to distinct points 

The aim of this section is to prove Theorem II. II in the case that all points 
are distinct and all bi as well. We first present the case r = 0. The multi- 
time joint gap probabilities will be discussed in the next section, implying 
the case of the one-time process; i.e., general r beyond r = 0. However, first 
presenting the one-time case will prove useful for understanding the basic 
structure. 

Denote by p{x, y; t) the one-particle Brownian motion transition from x 
to y during a time interval t, namely 

P{x,y\t) = -r=e 2t . (2.1) 

V iTlt 

Let us consider n + m Brownian bridges leaving at t = — 1 from flm+n < 
. . . < ttm+i < am < ■ ■ ■ < cLi and ending at t = 1 at positions hm+n < ■ ■ ■ < 
bm+i < bm < . . . < bi. The positions of these particles at time t are denoted 
by x(t) = {xi{t), . . . , Xm+nii)}- Then, the probability density that x(t) = x, 
conditioned that the Brownian bridges do not intersect in t G (—1,1), is 
given by Karlin-McGregor formula [30], namely 

P(x(t) = x) = ^ det{p{ai, Xj; 1 + t))i<jj<m+n det(p(xj, bj;l- t))i<i,j<m+n, 

(2.2) 

with Z the normalization constant, which is equal to the probability that the 
m + n paths do not intersect, given the initial and final conditions at t = ±1. 

It is known that a measure on x = (xi, . . . , Xm+n) of the form (12.21) has 
determinantal correlation functions (see e.g. Proposition 2.2 of [11], or for 
information on determinantal processes [7,25,31,40,41]). 

As mentioned before we restrict the discussion in this section to the case 
t = 0. Then, the /c-point correlation functions p^'^^ are given by 

p^^\xi,...,Xk) = det(ir(x„Xj))^<.^^.<^, (2.3) 

with the kernel K explicitly given by 

m+n 

K{x,y) = Y,Pi^^bf,l)[B-%p{a„y;l), (2.4) 
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where 

B = [Bij]i<ij<m+n, Bij= dxp{ai,x;l)p{x,bj;l). (2.5) 

Jr 

In particular, the gap probabihty of a set E, i.e., the probabihty that none 
of the Xi, . . . ,Xn+m belongs to the set E, is given in terms of a Fredholm 
determinant, 

P(none of the Xi G ^) = det(l -Xsi^Xi,)L2(R), Xsi^) = ^ E). (2.6) 

The structure of the measure does not change when taking the limit of 
one of more of the Brownian bridges starting and/or leaving from the same 
position. Thus the determinantal structure of correlation still holds, yielding: 

Proposition 2.1. Consider a^+i = . . . = a^+n = and &m+i = . . . = 
bm+n = and the other m Brownian bridges from ai to bi, with < < 
. . . < ai and < 6^ < . . . < 6i . Then, 

F{x{0) ^E) = det (1 - XEKn,mXE) , (2.7) 

where the kernel Kn,m is given by 

m 

K^,U^,y) = K^-^''''{x,y) + ^ 4-\x){f^-'%vf (y). (2.8) 

The Hermite kernel i^Hermite defined by the classical Hermite polynomials 
and their -norms^ 

71—1 ^ 

irr"'*'(x,i/) = e-(-'+^')/2 5^-/7,(x)if,(y); (2.9) 

i=0 

the functions ip^jf* and are defined as follows for 1 < k < m, 

„-x^/2 r -z'^+2xz p-xP- I'i C -z'^+2xz 



(2.10) 

where TQ a/2 denotes any contour containing the points z = 0, ai/2, . . . , 0^/2, 
and similarly for ro,b/2- Finally, the entries of the matrix of inner products 

/i = (/iM)i<M<m with fiM = {ip'k\^p^"'') = / dxi^'i^\x)tlj'['\x) (2.11) 
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can 



be written^ 



f^ke 



27ri 



dz- 



Z" [Z 



akh/2)' 



(2.12) 



Proof. We start from the setting ( I2.2l) -( l2l6l) and take the hmit when the 2n 
points ttm+n, • • • , o-m+i — * and bm+n, • • • , bm+i — > 0, and leaving the 2m 
points a^n < ■ ■ ■ < ai and b^ < ■ ■ ■ < bi fixed. Then, the probabihty density 
on the x,'s becomes 



P(x(0) = x) = -det 



v 



<j<m 
l<j<m+n 

l<j<m+n 



\ 



det 



)l<i<m » 
l<j <m+?i 



v 



l<j<m+n 



(2.13) 

where Z' is a normahzation constant. Consider any set of functions 
{lp'^^\x), /c = 1, . . . , n + m} spanning the vector space 



V{ai, . . . , Om) = span{e"*'' 1 < i < 



m, x^-^e~^ 1 < j < ^^}, (2.14) 



and similarly a set of functions {il)^^\x),k 
V{bi, ...,bm). Then, 



n + m} spanning 



^(x(0) = x) = i^det (y^f^l 



l<t,j<7i+m 



det{4-\xj) 



l<i,j<n+m 



[2.15) 



As mentioned above, this measure defines a determinantal point process with 
defining kernel 



n+m 



K{x,y) = Y,i^t\^)[B-\,vt\yl 



(2.16) 



where B = [-Bij]i<ij<n+m has entries Bij = {ipf'\ip''P^). Thus the goal is 

to find nice functions ■j/'^"'' and y?^"^ such that the inverse of the matrix B is 
manageable; usually one looks for a set of functions such that B becomes the 
identity matrix (bi-orthogonalization). In this instance, it is more convenient 
for doing asymptotics to find functions such that the matrix B has the form 



B 



/i 
1„ 



(2.17) 



^Similarly a^b^ denotes a contour containing and akhi/2. Note that 



dz 



2-Ki JFo.u z'^{z—u) u' 
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As will be shown below, the choice of functions for which this is the case is 
as follows: 

cpi' (x) = (p dz ; -— , 1 < k < m, 



An) 



(2.18) 

The Hj^{x) are the classical Hermite polynomials, with generating function 



j! ' 
(2.19) 



-^+2-- = 4^7(3;), and thus — / e-"'+2-" 

and with orthogonality relations 

/ dxHk{x)He{x)e-''^ = 5k,icl, with = ^2^^. (2.20) 
By the residue theorem it follows that 

d^^r ^ e span(e'^'=^ i^o, • • • , H^-i)e-^"/'. (2.21) 



Jro,a/2 ^"(^-afc/2) 



Similarly one defines the functions ip^^\x) upon replacing by 6^ in (12.181) . 
Thus the set of functions {(p^^\x), k = 1, . . . ,n + m} spans the vector space 
V{ai, . . . ,am), and the set {■?/'[."■' (x), A; = l,...,n + m} the vector space 
V{bi, . . . ,bm), as defined in (12.141) . 

The last step is to show that with our choice we actually obtain (12.171) . 
From the representation (12.181) of the (y9^"^(x), ^/'^"^(?/) in terms of Hermite 
polynomials, it follows immediately that 

fi^^ = T/^f = 6ke for m + 1 < k,e <m + n. (2.22) 

Next we show that 

{(p^^\ip^^^) = for l<k<m,m + l<i<m + n 

and m + 1 < k < m + n,l < i < m. 



16 



Indeed for 1 < k < m and ■m + l<£<m + n, we have 



'2 /. _^„2 



const / dz £ dw ^zw 



(27ri)2 Jp^^^^^ - afc/2) /r^ w 
const^vr / dz Pi_rn-l{z) 



■r , / X 0, (2.23) 

where Pi(a;) is a polynomial of degree i. The result is zero because for 
I — m — \ <n — \ the residue at infinity is zero. 

Finally, for 1 < /c, £ < m, by the same argument one gets 

(n) An)-, a/tt /" dz f dw e^"""" 



(2.24) 

By the residue theorem, the contribution of the pole at w = is a polynomial 
of degree n — 1 in z. Thus the integral over z is zero, because the residue at 
infinity is zero. Thus, it remains to compute the contribution of the pole at 
w = be/ 2, namely 

_ /• e'''' 2" _ 0F2" I ^ 

~ (2^^ 7r^_^^^^ 'z-{z-a,/2)^ ~ (2^ /ro,.,.,/. ^"(^ - a,6,/2) ' 

(2.25) 

This ends the proof of Proposition 12. 1[ □ 

The next step in showing Theorem 11.11 for r = is to determine the 
n — i> oo limit of the kernel under the space scaling 

with Oj, bi scaled as in fll.2p . 

a. = v^(l + ^) and 6. = (^1 - j (2.27) 
and with the assumption 

di <bj, 1 < i,j < m. (2.28) 
Thus, we have to show that for ^1,^2 in a bounded set, 

hm -^—^K^^^^(x,y) = ir^^(0;6,6)- (2-29) 



n— >oo 
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It is well known that the Hermite kernel under the above scaling, for ^i, .^2 
in a bounded set, converges to the Airy kernel Kj\^ (see e.g. App. A. 7 of [21]) 

=: ^^(^1,6), (2.30) 

where the path r> goes from e~^'^*/^oo to e^'^^^^oo, the path r< from 
to e~'^^^^oo, with r> and r< not intersecting each other. 

What remains is to compute the limit of the last term in (12. 8p . Since 
m remains finite, one can take the n —>■ 00 limit inside the sum. Below we 
compute the asymptotics for and n~j separately. Let us start with 

the matrix fi, as defined in (12.171) . 

Lemma 2.2. The following asymptotics holds for the inverse of the m x m 
matrix: 

lim J' , (—\ = -A-\ where A = ( — . (2.31) 



n—>co 



Proof. Using the scaling (11.21) . the quantity 



^ = " l + ^ + (2.32) 



^1/3 \^n^/3 

is, for n large enough, strictly less than n by assumption (I2.28p . We use 
(I2.12P and make the change of variable z = un 

du 



where 
with 



F{u) := u - \nu = 1 + ^{u - if + O {{u - if) , (2.34) 



Re(F(n)) = Re(?i) - ln(|u|). (2.35) 

Thus, we can deform the path \u\ = 1 into 75 = {1 + iy, —S<y<6} 
plus a circle segment 7' centered at zero joining the extremities of 75. By 
(I2.35p . the path 75 V 7' is a steep descent path for F with maximum at 
u = 1, F{1) = 1. We choose 5 = n~^/^, then, the contribution of the integral 
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in fl2.33p from 7' is of order 0{e ^"■^'''') smaller than the main contribution, 
coming from 75, for some c > 0. Thus, continuing fl2.33l) . 

IJ'ke= { — - — / du J — l + 0(e 

(2.36) 

By the change of variable u = {u — the last integral becomes 



du . . (2.37) 

27riy_,„i/io - (Ofc - 6^)^1/6 + (9(^-1/6) 



In the n ^ 00 limit we finally have 



lim V2ni/6(I237D = ^. (2.38) 



lim V2n^'^ - = r = (2-39) 

n^oo \2e/ ai, — ft* 



Thus, we have shown that 

n \" —1 

dk - bi 

This suffices to prove Lemma 12. 2^ since the dimension of the matrix does not 
depend on n. □ 

The next item is to determine the asymptotics of (p^^^ and ipk^^- 

Lemma 2.3. Consider the scaling 1^2. 26\) and ^HSj, with ^1,^2 in 0, hounded 
set. Then, 

(2.40) 

where Ta,,> is a simple path from e~^'^*/^oo to e^'^^/'^oo and passing onto the 
right ofdk- Similarly, 

MCl) ■■= lim (-) V-l"^ ( = — / du- ^, 

(2.41) 

where F^^^ is a simple path from to e '^^^^00 and passing onto the left 

of bk (similar to Figure\^ . 

Proof. The plan is to compute the large n behavior of 



ifl'Hx) = (p dz— — , (2.42) 
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Figure 8: Lines with Re{F{x + iy)) = Re(F(l)). 



with 



ak 



2n 1 



1/3 



n 



X 



2n + 



(2.43) 



Rescahng the integration variable z = u^/njl^ one gets 



2\ "/2 g-n-5ni/3+C'(n-i/3) 

n j 2'ni 



du- 



r , u — 1 — auln^/'^ 



(2.44) 



where F{u) = + 2u — In('u). The leading contribution comes from the 

neighborhood of the double critical point of F{u) at m = 1, where we have 



F{u) = \-\{u-lf + 0{{u-lf). 



(2.45) 



As integration path one can choose any path passing through u = 1 + an~3, 
with dk < ce, locally following the directions e^^'^*/^, and which remain inside 
the region G of Figure [HI Then, the integration away from a (5-neighborhood 
of M = 1 + an~3 (where 6 = , with < e < 1/3) will be of order 0{e'^^) 
smaller than the leading term, with < c ~ 5'^ for small 5. Then, in a 
5-neighborhood of m = 1 we can use series expansion and after the change of 
variable oj = n^^^iu — 1), one finds 



(n) 



2n 



/2e\ -/2 
\ n 



2T\i 



O 



n 



1/3 



duo- 



io - ak 



(1 + 0(e-^")) 
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where the integral goes from e '^'^^/'^Sn^/^ to e^'^^^'^dn^^'^ ^ and passing on the 
right of Ofc. From this, the n ^ oo hmit in fl2.4Up holds. 

The asymptotic for ip'i^\x) is essentially the same, except that uj h- > —uj 
and Sfc 1-^ —hki ending the proof of Lemma F2.3[ □ 

We shall also need: 

Lemma 2.4. Given the matrix 



A:=\\ 1 I , (2.46) 

a,- — 0,- 



the following identity holds: 



[A _ 1 ( Yj I w - ak 

Proof. Since 



y . =^ TT^^ -1 • (2.47) 



, , A(a)A(6) 

det A = ^ ' ^ \ , (2.48) 

one checks the identity fl2.47p . by computing the residue on the right hand 
side at the points z = di, w = bj and identifying with {A^~^)ij using Cramer's 
rule and repeatedly using (12.481) . □ 

Proof of Theorem \l.l\ Assembling the asymptotic result fl2.30l) . Proposi- 
tion Lemma YI7I\ and Lemma [2. 31 one obtains Theorem 11.11 in the special 
case r = 0, with distinct a^, foj, under the condition < hj. Upon using the 
scaling ( 12.261) and f l2.27p . the limit kernel is thus given by the limit of the 
sum of the kernels in (12.80 : i.e., the sum of the Airy kernel K^, defined in 
(I2.30p . and a new kernel: 

^ m 

= ^ / duj / du ^3 ^ (2.49) 



The fact that this expression actually equals the kernel i^'^''(0; .^i, ^2), as 
defined in Theorem ll.il follows from Lemma [2.41 □ 



21 



3 Extended kernel for the Airy process with 
wanderers 

In this section, we will prove Theorem 11.21 For this purpose, we need to know 
the measure, defined on the positions of the Brownian bridges at different 
times -1 < Ti < T2 < . . . < Ti < 1. Set x(Ti) := {xi(Ti), . . . ,Xm+n(Ti)). 
Then, by Karlin-McGregor applied to these different times, the measure ob- 
tained by the non-intersecting condition on the Brownian bridges is given 
by 

P(x(Ti) = x\ . . . , x(T^) = x^) = ^ det(p(ai, x], Ti + l))i<Lj<n+m 



X 



]^det(p(xf,X^+\Tfc+i - Tk))l<i,j<n+m I det{p{Xi,bj,l-Ti))i<ij<n+m- 
k=l J 

(3.1) 

It is well known that this measure, a generalization of (12.21) to multi-times, or 
any measure of this form has determinantal correlations in space-time [16,20, 
24,32,37] (even in cases when the size of the determinant is increasing [12,13]). 

Proposition 3.1. Any measure on {x'f'\ 1 <i < N,l < n < i} of the for\ 
det(0(To, ai] Ti,xf)i<i^j<N ( JJ det(0(T„, xf'] T„+i, xf''^^))i<i^j<N 

\n=l 

X det(0(T£, xf] T^+i, fej))i<ij<7v, (3.2) 

has, assuming Z 7^ 0, the following k-point correlation functions for 
ti, . . . , tfc G {Ti, . . . , Ti}: 

p^''\ti,xi, . . . ,tk,Xk) = det {K{ti,Xi]tj,Xj))^^^ .^f^, (3.3) 

where the space-time kernel (often called extended kernel) is given by 

K{ti,Xi;t2,X2) = -(l){ti,xi;t2,X2)l{t2 > ti) (3.4) 

TV 

+ ^ 0(ti, xi] T^+i, hi) [B~\j(l){To, aj; ta, X2) 

with (* means integration with regard to the consecutive dots) 

,(rj. rj. f 4>{Tr,x]Tr+i, ■)*■■■* (l){Ts^i,-;T„y), ifT^KTs, 

n-Lr,X,ls,y) - I ^jr^^ y rj,^^ 

(3.5) 



^The functions (j)(Tn, x;Tn+i,y) themselves may in fact vary with n above. 
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and with the NxN matrix B having entries Bij = (j){To, af, T^+i, bj) . Remark 
that [Niy det{B) = Z, so that B^^ exists as soon as Z ^0. 

We now apply this general fact to the non-intersecting Brownian motion 
formula fl3.ll) : here x-"^ denotes the position Xi{T„) of the ith Brownian 
motion at time T„, while one sets Tq = —1, T^+i = 1, and one sets 



(pit, x; t' , x) := p{x, x , t' — t). 



(3.6) 



As for the one-time case, the structure is unchanged, even after letting 



and 



0, keeping < . . . < ai and 



km < ■ ■ ■ < bi fixed. The only difference is that the entries on the first 
and last determinants in fl3.2p will be different (together with a different 
normalization constant Z). Indeed, the first determinant in (13.21) is just 
replaced by 



/ (e-4^V(i+TOp(o,xf,Ti + l))^ 



det 



V 



(tt7t)"V(o,x«,Ti + i; 



while the last determinant is replaced by 



/ feM''/(i-T.)p(^W^o,l-T,: 



det 



v 



(i) 

(^)-V(a:f,0,l-T, 



<i<m 
l<j<m+n 



l<i<n 
l<j<m+n 



l<i<m 
l<j<m+n 



(3.7) 



l<i<n 
l<j<m+n 



J 



As for the one-time situation, one looks for sets of functions generating the 
same vector spaces as the functions in (13.71) and (13.81) . namely one searches 
for functions (p^^\Ti,x) and ilj^^\Ti,x), such that 

dSH]) = const X det{ipf'\Ti,x]))i<ij<n+m, 
dM]) = const X det{ilj'f'\Ti,, x'j))i<ij<n+m, 



(3.9) 



and such that the matrix B has the same form (I2.17P as before. Setting 



l + t' 



(3.10) 



one picks, for 1 < k < m, 

-z2/2a(t)2 



x] 



a{t) 2m 
a(-t) 27ri 



dz- 



dz- 



z^{z-ak/2) ' 



(3.11) 



(z - bk/2) 
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and for 1 < k < n, 



(47r)V^ 



(3.12) 

Remark that, using the integral representation of the Hermite polynomials, 
an equivalent expression for fl3.12l) is 



(47r)^/'' e""" Z^''^-*) (A; - 1)! / g-^S(-t)'+2x^/'^2(_i) 



(3.13) 

It is immediate to verify that these functions generate at t = Ti, resp. 
t = Tf, the same space as the function in fl3.7l) . resp. (13.81) . So, one defines 
the functions appearing in the first and last determinant of (13.21) by 

0(ro,a,;Ti,x«) := ^(ri, x«) and 0(T,, x^; T,+i, 6,) := 4''\t,,x^'^), 

(3.14) 

for which we show the following property: 

Lemma 3.2. For any ti < t2 and 1 < k < n + m, one has 

dx ^f\ti,x) p{x,y]t2 - ti) = y?S^"^(t2,2/), 

(3.15) 

dyp{x,y\t2-ti)^^^\t2,y) = ip'j^\ti,x). 

Proof. Since ip^^ is obtained from y^^"^ by the map t —t and a i— 6, it 
suffices to present the proof for v^^^^ At first, for 1 < /c < m, one has 

dx(p''^\ti,x)p{x,y]t2-ti) = — (p dz- 



/• p-xV2(l+ti) -(x-j/)2/2{t2-ti) 

u v/mr V27r(t2 - h) 

and, after performing the Gaussian integration, one has 
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Secondly, consider 1 < k < n. Comparing the representations fl3.1ip and 
fl3.13p . we see immediately that the computations are exactly the same. In- 
deed, the only difference is a fc-dependent prefactor and the denominator in 
the integrand over z. However, there are not affected by the computations 
above; thus 

/ dx^';^l^{ti,x)p{x,y;t2-ti) = ^l^li^{t2,y) (3.18) 

holds, ending the proof of Lemma 13.21 □ 
Proposition 3.3. The extended kernel is given by 

n 

Kn,m{tl,Xi]t2,X2) = -p{Xl, X2; t2-ti)l{t2 > + Xl)Vm+iit2, X2) 

1=1 

m 

+ ^l,t\h.x,)[^Ji'\,vf\t2,x2). (3.19) 

with (p^J^\t,x) and ilj^^}^^{t,x) given by ^3.11\) and ^3.13\) and with n given by 
1^2.24^ , the same as in the 1-time case. 

Proof. Given the definitions (13. 6p and (13. 5p , the first term in the kernel (13. 4p 
is simply 

-(P{ti,Xi-t2,X2)l{t2 > ti) = -p{xi,X2,t2-ti)l{t2 > tl). (3.20) 

It remains to be shown that B has the form 

as in the 1-time case, with /i given in (I2.24p . 

Indeed, for any choice of 1 < /c < £ — 2, and for ti = Ti with 1 < i < i, one 
has, using the convolution property of the Brownian transition probability 
and the convolution property in Lemma 13.21 the property that (* means 



25 



integration with regard to the common variable) 

Bij = (f){To,ai;Ti^i,bj) 

= y.S"Hti,x«)*p(x«,x(2);t2-ti)*... 

...*p{x^'-'\x^')-t,-t,_,)*i,f\u,x^'^) 
= [y,t'\t,,x('^)*p{x('\x^'^'^;t,^,-t^)) 

*(p(x('^+^),xW;t,-4+i)*^f)(t,,xW)) 
= ^fHt,+„x(^+i)) * V^f (4+i,x('=+i)) = (v^S-H^+i, .)) 

is independent of tk+i] therefore, by setting tj^+i = 0, it is, in particular, 
equal to the value fiij obtained in f|2.25p and (12.241) . This establishes Propo- 
sition 13.31 □ 



In order to prove Theorem 11.21 (and thus also Theorem 1 1.1 1 for generic r), 
one needs to compute the n — > oo asymptotics of the kernel. For convenience, 
recall the scaling for the starting and ending points of the top m Brownian 
bridges (11.21) and of the subsequent scaling (II. 3p of the space-time region one 
focuses on: 



'2n + 



,n ii-^f (3.22) 

72^1/6 

with di < bj, 1 < i,j < m. Below we prove that, given the scaling (13.221) and 
for ^1,^2 in a bounded set, 

hm ^^i^„ ,„(ti,xi;t2,a;2) = K^^(ri, 6; r2, 6)- (3.23) 

where = we means an equivalent kernel^. 

Proposition 3.4. With the above scaling, for ^i,C,2 in cl bounded set (and 
Ti,T2 fixed), in the case where all the di (and bi) are distinct, one has 

hm ^^ir„_^(ti,a;i;t2,a;2) = i^^'(ri,ei;r2,6)7r^4' (^.24) 

where /(r, ^) = exp(r^/3 — ^r). 

^Two kernels are equivalent if they define the same determinantal point process. 
Namely, if there exists some function f{x) ^ such that K{x,y) = K{x,y)f{x)/f{y), 
then K and K are equivalent, since all the correlation functions are given by determinants 
in which the functions / cancel exactly. 
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Proof. Consider the first two terms in the kernel (13. 191) . These terms are 
independent of the a,, bi and of m. Indeed, it corresponds exactly to the kernel 
of the system without wanderers, which can be denoted by Knfl- Indeed, 



k=0 



(3.25) 



For fixed ti,T2, we show below that 



uniformly for ^1, ^2 in a bounded set, with the extended Airy kernel given 
by 



KA{ri,^l]T2,^2 



L dX e^(-^--i) Ai(ei + A) Ai(e2 + A), n > r2. 



- 4_ c/A e^(-2--i) Ai(ei + A) Ai(6 + A), n< r^. 

(3.27) 

To obtain this result, for ^1, ^2 in a bounded set, one can just use the asymp- 
totics of the classical Hermite polynomials (see for example App. 7 of [21]). 
Another, better, way is to first perform the sum over k using two different 
integral representations for Hermite polynomials, a first one is fl3.13p and a 
second one is an integral over L + (see e.g. sect. 2.2 of [29]) for L > 0; 
namely: 

Then 

Kn,o{Xi,ti;X2,t2) = -p{Xi,X2;t2 -ti)l(t2 > ^l) 

^1 ^2 /TT\n l-ti rr2 2xiU 

+ ..2 =^ / ^ ^^ ^r/ ^37--^. (3.29) 



(27rz)V(l + ^i)(l + ^2) A+« A U-V ^^v^-'-Bf, 



Note that the — 1 in (^)" — 1 (appearing in the integral above) can actually 
be omitted, because there is no residue aX V = U . One then makes the 
substitution to new integration variables f/ and V , 
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and uses steepest descent in the integral to get the extended Airy kernel 
fl3.27p . which is just (ll.7p in which one replaces dk = bk = and m = 0, 
namely 

KAru(ur,,(,) = - '("^ > -g^fai;-|(^-n)(6.6HA(-,-n,3 

V47r(r2 - Ti) 

+ - -5 dto duj — ~3 in^c ~ 1 ^ 7^ 7- (3.31) 

What remains is to compute the limit of the third term in fl3.19p , namely 
1™ -T^e E 4''\h,x,)[f,~%^f{h,x,). (3.32) 

n->oo -i/2nV6 

Since m remains finite, we can take the n ^ 00 limit inside the sum. Also, the 
limit of taking into account the prefactor, has already been computed 
in Lemma [2.21 It remains to determine the asymptotics of ilj^\ti,Xi) and 
ip^j^\t2, X2) (for 1 < k < m) under the above scaling. 

As will be seen, the computations are very close to the ones for t = in 
Lemma 12.31 For convenience, recall the notations 7(t) = a/(1 — t)/(l + t) 
and a(t) = \/l + 1. From (13.111) . after the change of variable z = w/-yit), 
one gets 



where x' and are defined below, together with their asymptotics 

X X r— ^ 



a'^ ■= akiit) = + V2{dk - r)n^/^ + 0{n-^^^). (3.34) 



Now we benefit from the computation made in the r = case. Indeed, we 
showed that 

^^7 ^ = /M — — / du ^ (1 + 0(1)) , 



2'Ki An , z"'(z — a/2) \ n / 27ri Jr- 00 — a 

(3.35) 

if y and a are scaled as 

y = V2n + ^ a = V2n + V2dn^/^. (3.36) 

v2 n^/^ 
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This is exactly our situation with a = aj^ — t. Thus we get 



ip^^\t,x)= 7(t)"e"/M- — / du- (l + o(l)). 

(3.37) 

Moreover, the asymptotics of the prefactor reads 



ait) 

Thus we have showed that 

n— >oo \ 

and similarly, 



„-x2/2(T(t)2 



doj- 



2e/ 27ri 7ra,_.2> ^ - ^fc + ^2 /(r2, 6) ' 

(3.39) 

(n) (n\-^l'^ 1 /" ~ (p'^l'^-^^'^ 
^kijx.ii) := \\m ^]^\tx,xx)\ — ) =— / f /(n,6)- 

n^oo Vie/ 2712 _ ^ _ 5, _|_ T 

<''fc-n 

(3.40) 

Now we can put together all the pieces, which make up the kernel (13.191) . 
namely (I3.26p . (13.391) . (I3.4UI) . and the asymptotics of the inverse of the matrix 
B in Lemma [2.31 Thus we have 

- E Mri,Ci)[A-\mir2,C2). (3.41) 

The last term in (I3.4ip (including the minus sign) is equal to 



{2my Jr,_.,> Jr,,_^^ e-^/^+«^- _ 5, + n) (^ - a, + r2) 

(3.42) 

Applying the identity in Lemma 12.41 we get as final result the kernel 
-^m''(''"i5 ^1! "^25 ^2) of Theorem 11.21 and this ends the proof of Proposi- 
tion El □ 

Proof of Theorem \1.2: For any bounded set E, the probability fll.lOp is given 
by the Fredholm determinant of the kernel, obtained in Proposition 13. 4[ 
Since this kernel is conjugate to the one in Theorem 11.21 their Fredholm 
determinants are identical. □ 
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Proof of Theorem \1.3\ (Universality). The proof is a mild variation on the 
proof of Theorem 11.21 and Proposition 13.41 Referring to the notation used 
in the statement of Theorem 11.31 one checks that formula (12.321) for akbi/2 
remains the same, since XqXq = 2n (see (II. 8p ) and thus also asymptotic 
formula (12.331) for ^ki- Moreover the scaling now reads: 



X 



VMl-t') (l + ^) (3.43) 



tor toT^ 



2„(l-«(l-^ + ^-|^K + r') + 0(„-/3: 

Refering to the notation (I3.34p . one checks that with x and t as in (13.4; 
above, one has 

r' = - = - = ( 1 ^ ^ \ M nir,'^l^\ 



x' = , , , = , = V2n 1 + — ^ + Oin'^'n 

a'k = aklit) = ak^\^^ = (^1 - + 0(^-1/^) (3.44) 

With this information, one checks the following asymptotics, which is the 
analogue of (13.381) . namely 

^ -7(±t)"e^ /2 e 2 /i n 



a{±t) " ' VT±t\l + t 



with 



(3.45) 



Here fn{'T,C,) depends on n, besides r and ^. Then we show 



— / du——^ , (3.47) 

27rz Jr. uj- ak + T2 
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and similarly, 



/ du—^ . (3.48) 



— 1 /o \ " 



Also, as before, 

lim ( — ) = -A-\ (3.49) 

Then, with as in 03.431) . the limit (I3.26P gets replaced by 

^™ fK -^n,o(^i^^i;^2,a^2)^ — ^ = ^^(n,4i;T"2,42), 

with very little change in the steepest descent argument. So, putting all the 
pieces together, one checks: 



/F: 1/. I<-n,m{tl,Xut2,X2)—. = A (Ti , 4l ; , ^2) 

n^oo A/2nV6 Jn (Tl,4lj 



(3.50) 

from which one proceeds in the same way as in the proof of Proposition 13.41 
and Theorem II. 2[ This ends the proof of Theorem 11.31 



4 Airy process with wanderers all leaving 
from point a and all going to point b. 

In this section we prove Theorem 11.21 (and thus also Theorem 11.11) for the 
case where m wanderers all leave from one point and all are forced to one 
point; i.e., 

d := dm = ■ ■ ■ = (ii < bi = . . . = bm ='■ b. (4.1) 
Thus, the m top Brownian bridges start from a and end at b with 

a = V2^{l + dn-^/^), b = V2^{1 -bn-^^^). (4.2) 

The arguments presented in the previous sections break down. Therefore 
one should redo the proof, using an argument adapted to this case. It is 
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instructive to shortly present two different approaches. The first follows 
the approach of the previous section, consisting in computing the inverse 
of the m X m matrix /i, and the second approach is to perform the bi- 
orthogonalization. In principle, with some care because of the n ^ oo limit, 
one might also be able to do the argument by analytic continuation, since 
the measure is analytic in the well as the final kernel (provided the 

inequality di < bj for all i,j is satisfied). 



4.1 Via the inversion of the moment matrix 



The start is almost the same as in the previous section. The only difference 
is that the first and last determinant in the measure, instead of (13. 7p and 

), are now 



det 



^ ((x])'-v^.'/(i+^^V(o,x],ri + r 

((a;])^-yO,a;),Ti + l))i<.<„ 

l<j<m+n 



\ 



l<i<m 
l<j<m+n 



(4.3) 



and 



det 



(x^)*-ie'"'/(i-^^W,0,l-T,) 



v 



((xj)*-V(xj,0,l-T,))i<i<„ 



l<j<m+n 



l<i<m 
l<j<m,+n 



respectively. The functions v^^"'' and for 1 < k < m, defined by 



(4.4) 



a{t) 2m 
a{-t) 271^ 



dz 



^-z^-y{tf+2xz/(7^{t) 



0,a/2 



dz- 



0,6/2 



{z - a/2) 
(z - hl2f 



(4.5) 



replace those of fl3.1ip . where we recall that 7(t) = y jq^' ^(^) ~ + 1. 
Of course, since the last n rows of the determinants (14.31) and (14. 4p are exactly 
the same as in (13. 7p and (13.80 . we keep the same choice for the functions v^^^^ 
and V'm-i-fc' 1 < < as in (I3.12p - (l3.13p . Define the matrix m x m matrix fi 



by fiij = {ip^P\ilj^"^) , I < i, j < m. Once again, this choice of (p)^'' and tp^ 
generates the same vector space as the function in the above determinants 
gSD and diaD. 

Note that in this section we use the same notations as in the previous 
section. However, the matrix /i and some of the functions are not the same. 



(n) 
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What remains the same is the form of the kernel. Indeed, since Proposi- 
tion EiS] holds exactly as before, one has once again: 

m 
i,j=l 

(4.6) 

where is the kernel without wanderers and in the scaling limit will 

converge to the extended Airy kernel. Thus, we only have to deal with the 
double sum below. 

Lemma 4.1. Under the scaling ^4-2^ , we have 

/ \" / , 1 1 /l + k-2\ , , 



/or 1 < k,l < m. 

Proof. For convenience, in the proof we compute fik+i.i+i to avoid — I's in the 
formulas. Since, as before, fike is time-independent, we may set ti = ^2 = 
in the computation; so, as in (12.241) and after integrating over the x variable, 
one finds 



Then we apply twice the identity 

1 2'' f d\'' 1 , 

(4.9) 



(z-a/2)'=+i k\ \daj z - a/2 
and obtain 

^^^^''^^ = WUJ w ^'-'^^ 

But /ii^i was already expressed as a single contour integral; see (I2.12p . Thus 

= W (a; j (i; j ^ ''^ ."(.-aV2)' 

We now compute the derivatives of {z — ab/2)^^ and obtain 

2^ 7^ Y ^ ^ , "^^-'^ (a6/2)-^- (/ + fc-j)! 

fc!/! \daj \daj z-f ^ (^ - f )'+^-J+M'^ - j)K^ - ' 

(4.12) 
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and so 

min(A;,Z) 



x^^d> dz— , , (4.13) 

Finally, we need to do the asymptotic analysis of the integral, which 
essentially has already been made in Lemma Consider the following small 
change in fl2.33p . with = a, bi = b: replace {z — ab/2)~^ by {z — ab/2)~'P~^ 
for any finite p = 1, 2, . . .. Then the steepest descent analysis is unchanged 
except for that finite power, which would be present in (12.371) too. This extra 
power gives a factor n~^/^ and we also have an extra factor coming from the 
change of variables equal to n^*'/^. In the end, the result is 

7r2" P /2p\" 1 r)-2p/3 

d^—, ^1 =1-) — ^ ^^—(1 + 0(1)). 



J^oa„2 z^{z-ab/2y+^ \n ) (6 - S)p+i v^nVe 

(4.14) 

We put fl4.14p into (14.131) and compare the dependence in j of the different 
terms in the sum. Since ab/2 = n(l + 0{n^^^^)), the j'th term in the sum 
(I4.13P contains the following power of n, namely: 

Therefore, since the sum is finite, in the n oo limit, the leading term is 
the one with j = 0, the other ones being of smaller order. Thus, 

fk + l\ f2eY 1 a^fe'^n-^C^+O/s ^ 
1 fk + l\ f2e\ / V2 



, , — — (1 + 0(1)), (4.16) 

ending the proof of Lemma 14. 1[ □ 



Corollary 4.2. With the same scaling as in Lemma \4-1\ we have 

( \ k+l-l 
=2(i-ar'^-^HL-YL-%, (4.17) 

for 1 < k,l < m, where the mx m lower-triangular matrix has binomial 
entries 

(^-^).,/=(-i)'^-'(^r^^)- (4.18) 
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Proof. From fl4.7p it follows that computing the inverse of fi reduces to com- 
puting the inverse of the m x m matrix u, 

= ii^k,i)i<k,i<m, with z/fc^i = ( ^ _ ^ j . (4.19) 

A convenient way of taking the inverse is to compute the z/ = LL^ decom- 
position, where L is lower-triangular, yielding 

u = LL^, Lk,i = r " M , (4.20) 



l-l 



from which 



r^-' = {L^rL-\ iL-\,i = i-^r'[]_l)^ (4.21) 

This establishes the asymptotics (14.171) . □ 

We now turn to the asymptotics of the functions v^^"-' and defined 
in (ESI). 



Lemma 4.3. Under the scaling l[4.2^ and 

c. _ 0-2 

t. = nn-^/^, Xi = V2^ + — (4.22) 

x/2nV6' ^ ' 

one has 

n^oo \Ze/ 



rfc^e"" /'+^^" / ^ "'^"^ (4.23) 
27ii Jr- iu! — a + T2)'' 

^ a — T2> ^ 



and 



n^oo \Ze/ 
^ i-iy^ r ^^g^V3-6^ /(^I'^i) (4.24) 



<b-Ti 

Proof. We must compute the asymptotics of 



^pU,x) = - -6 dz- ^ (4.25) 
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and compare this expression with fl3.11l) . One sees that the only differences 
are that now a replaces and the denominator in 2; — a/2 has a power k 
instead of power 1. For any finite k, the asymptotic analysis for this case has 
only minor differences with respect to the asymptotic for (13.111) . Namely, one 
picks up some extra factors by the changes of variables: setting z = 

one gets a factor y2/n and then to = n^^^{u — 1) results in a factor 
^('=-i)/3_ total, an extra factor {\/2/n^^^)''~^ appears. A similar argument 

holds for V'l"^- n 
Proof of Theorems \l.l\ and for Brownian bridges starting from a and 
ending up at b. Putting together Corollary 14.21 and Lemma 14. 3[ one obtains 

^™ /TT 1/fi i^i''\tl^^l)[f^'\jVf'\t2,X2) 

n^^o V2 n'-'^ 

= E^.(n,ei)(&-s)^+^-M(^^)~^^-^]..^.(r2,6)7^ 



1 



/(r2,6) 



(27ri) 



Fa 



-T2> 

m m 



duj ■ 



-a;3/3+^2t^ 



/(n,ei 



/3+6^/(r2,6) 

(I _ a)^+J-2(_l)i 

V" J -ij (u + n- by{uj + T2- ay 

Finally, using the fraction decomposition identity 
1 //f/_a\" fV-b " 



(4.26) 



<h-Tx 

k-\\ (k-\ 
I - 1 



u -b 

m m 



1 



V -a^ 

k-i\ fk-i\{b- ay+^-^{-iy 



(4.27) 



{u -by{v ~ay 



k=l i,j=l 

with U = uj + Ti and V = u + T2, one gets the final result 



1 



(27rz) 



Fa- 



^2> 



{u - a + Ti){uj - b + T2) 
(cD - 6 + Ti) (cj - a + r2) 



/(ri,6) 
/(r2,6)' 



(4.28) 



4.2 Via bi-orthogonal functions 

Here we present a slightly different approach, which consists in using 
biorthogonal functions, instead of the functions defined in (14. 5p . We use 
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a representation with determinants known from classical orthogonal polyno- 
mial theory. Let us first define the polynomials <^^"\ ^i"^' ^^"^ then show 
they are actually biorthogonal. 
Set 

/Xjj(t) := (V5f''(0, ■),^j"^(0, ■)) and := det(/iij)i<ij<fc. 

and define 

/ /^l,2 • • 

1^2,1 1^2,2 ■ ■ 



(4.29) 



1 



XI 



det 



fc-i 



and 



fik,k-i ipP{t,x) ) 



(4.30) 



det 



VV^i"^(t,x) 4")(t,x) 



fJ'2,1 



in) 



(4.31) 

First of all, notice that (p^^^ is linear combination of the (p^f^ with i = 1, . . . ,k, 
with a non-zero coefficient in front of ip^^^ (because ^ 0, since both 
{(p^^\l < i < k} and < £ < k} form a basis of a /c- dimensional 

vector space). The argument is similar for "ipl"^ . Observe, for i < k, 



right hand side of fl4.30p with 
the last column replaced by 
the £th column of fOOl) 



0. (4.32) 



Therefore also {p^j^\t, ■)^il)\"''{t, ■)) = for £ < A; and thus also for £ 7^ A;, by 
merely interchanging the roles of (p and ip. The above argument also shows 



-,(") 



(4'^^(t,-),C^(t,-)) = (^r(^r),C^(t,-)) 



AfcAfc_i 



1. 



The consequence is that now the kernel instead of (14. 6 p reads 

m 

-,-,-11 r»i 



(4.33) 



Kn,ni{ti,Xi]t2,X2) = K^fiih, Xi, t2, X2) + ^ ^f"'' (ti , Xi ) [/i \jpf'{t2,X2) 

i,j=l 

(4.34) 
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with = {ifi' \t, \t, ■)) = 5ij. Therefore, the double sum in f l4.34p 
becomes just 

m 

Km{tuXi;t2,X2) = tpl'^\ti,Xi)[fl'%jlp'f\t2,X2) 

m 

= Y.i,t\h.x,)^t\t2.X2). (4.35) 



This last sum is of Darboux-type and can be rewritten as 

( 4"^(ti,xi) ■■■ ^t\t^.x^) \ 

<~P^f'{t2,X2) 



A. 



■det 



\^ i{>:m'{t2,X2) 



/^2,1 



/^1,2 
/^2,2 



/^m,2 



(4.36) 



the latter follows from the fact that Km{ti,xi;t2,X2) is a bilinear combi- 
nation of iljl^\ti,xi) and (y9j"''*(t2, 2:2), for 1 < i,j < m and is completely 

characterized by {Km{ti,xi;t2, .),il'i^\t2, •)) = '4'i^\ti,xi) for 1 < i < m. 

At this point we have to determine the n ^ 00 limit of the rescaled 
kernel, namely 

lim -^=——Km(ti,xi;t2,X2) (4.37) 

with Xi,ti scaled as (14.221) . The asymptotics of fiij already appears in 
Lemma 14.11 and for yj^"'' and tpl"^^ in Lemma 14.31 Hence Lemma 14.11 to- 
gether with the fact that m is finite, yields the asymptotics of A^: 



/ fl \ nm 

lim Am — 

n— >oo \/e. 



n 



1/6 



V2 



1 



2™ (b - a)"' 
1 1 



det 



k + l~2 
k-1 



2™ (b 



l<k<m 

(4.38) 



This result (14.381) and the linearity of the determinant, together with the 
results of Lemmas 14.11 and 14. 3[ substituted in fl4.36p lead to the limit in 
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fl4.37p . namely 



a-T2> 



duj I duj 



g-cj3/3+^2<^ 



-(6-5) 



X det 



1 

1 



{U + T2 - a) {u + Ti - h) 

m— 1 \ 



( b-a 
\uj+T2—d 



\ k-1 ) l<k,l<m 



a;+ri — 6 



m—l 



/(t2,6)' 



(4.39) 

The final step is to use the following identity, established by observing that 
both sides are identical upon integrating against 1 < i < m from x = 
to X = 1, 



det 



/ 
1 



1 



{x + l)™~i \ 



l<kl<m 



m—l 



{xyY 



(4.40) 



with y + 1 = {a — b) /{u + Ti — b) and x + 1 
one obtains for f l4.37p : 



Km{ti,Xi;t2,X2) 



xy — 1 
a) / {uj + T2 — a). Thus, 



/ 



lim —= — - 



duo 



duj ■ 



-T2> " <b-Ti 



-a;3/3+6'^ 



X 



(b-d) 



W + Tl-b LU + T2—d 



(UJ + T2- a) (U + Ti- b) a;+ri-a u^+Tj-b _ 



(4.41) 



which is equal to the kernel (14.281) obtained previously. 

Remark: Using fl4.30p . (14.310 . a 1-border identity analogous to the 2- 
border identity (14. 40 p . Lemmas 14 . 1 1 and 14 . 3 1 and (I4.38p . one finds the limiting 
biorthogonal functions (which also yield (I4.4ip ): 



lim 2^/V/i2^"(t,a;) = Vb^L^)^[ duoe-^^^^ 

lim2VVA2^^(t,x) = Vb^atDLl f rfa;e^'€- (^ + ^~^j''\ 
n^oo -J 2'ni {u + T - by 



^^+^Juj + r-b,y-^ 
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5 Limit to the Pearcey process 

In this section we prove Theorem 11.51 To do this, we must apply the scahng 
(11 .lip to the kernel with m wanderers, where all the di = d and hi = h and 
where one uses the shift uj ^ oj — T2 and u; — > — ri, to yield: 

i<-S;'(ri.f.;r..&) = - > ^1) -t^-i(.-n,(6.6).A(^-n)- 

A/47r(r2 - Ti) 



LU — a 



(5.1: 



According to Theorem 11.51 the scaling limit we need to take is the following: 

a = am^^^, b = jSm^^^, 

Ti = T^m^''' + ]^^eim-^'\ (5.2) 

ii = Xm^l'^ - K^a±e,m^/^ - KVim-^/^\ 

Then, we have to compute the large m limit of the rescaled kernel (15.11) . We 
prove the following result, which implies Theorem 11.51 

Proposition 5.1. Under the above scaling, for any fixed 61,62, the limit 

hm Km-'/''K^\T,,^,;r2,^2) = K''{6,,v,;62,V2) (5.3) 

m— >oo 

holds uniformly for Vi, V2 in a bounded set. 

Proof. The first term in (15. ip is a straightforward limit. Indeed, for 62 > 61, 



1 (^2-n)" Q(1) ^ ^ 

e 2(92-61)^1^ (5.4) 



'2n{62 - 6,) Q{2) 
where the conjugation terms Q{i) are given by 

Q{i) = exp + X)6,m'/^ + naVim^l^ + 0{m-^l^)^ . (5.5) 

Next, one deals with the double integral, where it is natural to introduce 
the change of integration variables: 

10 = wm}^^ , bj = wm}^^ , (5.6) 



40 



leading to 



dw / dw 



(5.7) 

where the functions Fi are given by 

Fo{w) = -hw - Tf + X{w -T) + ln{w - P) - Hw - a), 

F^iw, 9) = l{iw- Tf - X)K'e -{w- T^ae, 

Fs{w, v) — —{w — T)kv, 

F^iw, 6) = ~{w -T- a)K^e'^. 



Setting 
one defines 



w 



':=w-T, a:=a-T, p' = (3-T, (5.9) 



Fo{w') := Fo{w' + T) = -— + Xw' + log(w' - a') - log(w' - f3'). (5.10) 

One now imposes the condition that F'{w') experiences a triple zero at some 
critical point w'^; this happens when the following polynomial P{w') is iden- 
tically zero, with w'^^ w'^: 

= P{w') := -{w' - a'){w' - /3')^oK) - W - OH^' - w[) 

= (3^; + w[ -a'- (3')w'^ + {a' (3' -X- 3wf - 3w>;)w'^ 

+ {w'^ + ^w'^w[ + X{a' + P'))w' - {a'p'X - a' + p' + w[^w\). 

(5.11) 

Setting the coefficients of this cubic in w' equal to amounts to 4 equations 
in 5 unknowns a', j3', w'^, w[, X, thus yielding an algebraic curve. At a first 
stage, let us look at it purely algebraically; later we will have to take into 
account the real character of the parameters, including various inequalities. 
Close inspection of the four equations suggests the following birational map 

""-^23^' ("2) 

with inverse (assuming a' + /?' 7^ 0) 

1/^/ A 2{2w[ + a' + f3') 1 , , 

r-^{P-a), 3^^^^^^^ , a^--{w,-a-P). (5.13) 
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Substituting this map in P{w'), one now solves the 4 equations (15.111) defined 
by P{w') inductively, beginning with the highest degree in w'. At first, one 
checks P{w') = 3{w'^ — a)w'^ + leading to w'^ = a, together with the value 
of w'l, which we already knew from (15.121) : thus 

w'= a and w\ = — -. (5.14) 

2 - X 

Substituting this back into P{w') yields a quadratic polynomial in w'] the 
vanishing of the coefficient of w'"^ yields 



2(T2(3a;2 _ 6x + 2 

(x-2^2 



X = S 7^, --^-^ (5-15) 



P{w') becomes thus linear, with vanishing linear and constant terms, yielding 

2 2 2 (23^-3) 2a;V 
^ =^ ^ (a; -2)2 '''''^ ''"2^- ^ ^ 

The compatibility between the r and equations (I5.16P yields a curve re- 
lating X and 0", 

(5.17) 

Incidentally, this curve is elliptic; indeed, viewed as a 6-fold cover of the 
x-plane, the total ramification index equals 12, with a ramification of index 
5 above x = 3/2, there are two ramification points of index 2 above a; = 
and three simple branch points above x = oo; thus the genus = 1. Then 
substituting the value (15.161) of r into (I5.12p and (I5.15p . yields the following 
expressions for a', [3' and X, all defined on the algebraic curve (I5.17p : 

a' = ^^(l-xV^') , /?' = ^^(l+a;V), and X = ^'(1 - 2a;). (5.18) 
2 — X 2 — X 

Using these expressions, together with the value of the critical point w'^ = a, 
one checks from (I5.10p that 

and thus 

Foiw') = Foiw',) + ^-^^iw' - w'y + O {{w' - w'^f) . (5.20) 

One then requires the parameters a', X, w'^, w[ to be real with a' < (3' , 
w'^ 7^ w'l, and a' + {3' 0. This implies that x, a and r must be real; the 
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curve relation fl5.17p yields two real solutions for a, namely o"+ < and 
cr_ = — cr+ > 0. In particular from fl5.17l) one must have x > 3/2, and 
since a' < one must have, according to (15.131) . that 2r = /?' — a' > 
yielding from f l5.16p the inequality x < 2. Thus one has 2 > x > 3/2. 
Moreover f l5.19p will be < for a = a+ and > for a = a-; one then 
sets =F^ := ^r^, with k > 0. Since ^f^, the right hand side of (15.171) . 
is an increasing function of3/2<x<2, this function has its maximum at 
X = 2, for which a± = =Fl/2, according to the curve relation (15.170 . That 
a' + (3' ^ follows from adding the two first equations in (15.120 . Also one 
has w'^ = a+ > [3' and w'^ = a_ < a', since from (15.181) and the curve (15.170 . 
one computes for (3' — cr+ and similarly for a' — cr_ : 



/3' - a+ = ^^(1 + 2xVi) = - ^2x - 3) < and a' - a_ > 0. 

2 — X 2 — X 

Finally, from (I5.14p . it is clear that w[ ^ w'^, since °^^2-x^"* 7^ in the admis- 
sible range x G (|, 2), with w[ — >• w'^ for x — >• |. 

To summarize, using the change of variables (15.90 . the relations (I5.18P 
imply, for a given a < (3, two values T± of T below, and thus 



« = T± + |^(l-a:Vi), /3 = T± + |^(l + xVi) 
X = 4(l-2x), T,= " + ^ 



(5.2i: 



2 2-a;' 

from which (I1.14p in the statement of Theorem 11.51 follows, with inequalities 

3/2<x<2, 0<|cr±|<^, X < 0, T_ < < T+. (5.22) 

Also the critical point Wc of Fo{w) and the extra-root wi of Fq{w) occur at 

, a<(3<Wc a±{3x-2) 

Wc-=(T± + i±, With and Wi = I±-\ Wc- (5.23) 

Wc < a < p 2 — X 

The statement about the uniqueness of the solution {x, a) to the equations 
(I5.2ip and (I5.17p . given arbitrary a < (3, remains to be shown. Indeed, upon 
using the identities (15.210 obtained for a and P, together with the curve 
equation (I5.17p . it is easy to see that the right hand side of the equation, 

0<(3-a = -^ = )- 5.24 

2 — X 2 — X 

is a monotonically increasing function in the range 3/2 < x < 2; therefore 
the right hand side of that equation takes on every value in (0, oo) exactly 
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once and thus given arbitrary a < /3, there is a unique value x G (3/2,2) 
satisfying the first equation in fl5.24p . Substituting this value of x in the 
T±-equation of f l5.2ip . the value of T± is specified unambiguously and thus 
T± can take on any value in M. Therefore, only when a, /3 — 0, do 

X 3/2, cr± — s> and T± — >• 0, and thus, by (15.231) . Wc^wi — > 0, (5.25) 

which proves the remark at the end of Theorem 11.51 

Then, the series expansions about the critical point a± give 

F^iw, e) = F2{w„ 6) + ]^K''e{w - We)^ (5.26) 

F^{w, v) = Fs{Wc, v) - Kv{w- Wc), 
F^{w,0) = 0{w-Wc). 

We now apply the steepest descent method, which we spell out for the open- 
ing cusp; i.e., for T+ and a = 0"+ < 0. By Cauchy's Residue theorem, one 
can deform the paths as indicated in Figure M The contribution of the last 
contour is zero. Indeed, the integration over w is trivial, since the only pole 
is simple at w = w. All the factors involving a and f3 cancel exactly. Thus, 
we remain with a contour in w around f3 of an analytic function (no pole at /3 
anymore) which is zero. The deformation also involves contributions which 
vanish at infinity. 

The final and most important step is to deal with the previous last con- 
tours of Figure First a remark on the integration paths in (15.71) . For 
large w and w, the leading term is the cubic in Fq, which means that with- 
out any error, we can let the the directions of the path w go to infinity in 
the cones of angles in {tt/2,5tt/6) and (— 57r/6, — 7r/2) instead of 27r/3 and 
— 27r/3. Similarly for w we can let it go to infinity in the cones with angles 
in (7r/6,7r/2) and (— 7r/2, — vr/6) instead of 7r/3 and —ir/S. Finally, the small 
contour around /? can be also deformed to go to infinity as soon as it does 
in directions in (57r/6, 77r/6). Therefore without errors we can deform the 
contours to become as in Figure [THl Let us verify that these paths satisfy 
the steepest descent property. This will be done for the case a = a^; the 
case cT = cr_ is essentially the same. 

Slope of the function Fq{w) starting from Wc = (t^. Consider the curve given 
by 

w = Wc + (^^e±^("/2+^), for < (5 < 7r/3 and C > 0. (5.27) 
2 — X 
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Figure 9: First deformation of the paths which then pass close to the critical 
point Wc = cr_,_ + T. The solid contours are for w, while the dashed ones 
for w. (For the case of (7_, one has Wc < a and the figures is essentially 
reflected.) The circles in the flrst and second flgures on the right hand side 
become dashed and in the second flgure it sits on the right hand side of the 
full curve. In the third figure, the inner circle is dashed and the outer is full. 
In effect, the roles of w and w' are interchanged. 



Remember that u < 0. Then, at flrst one verifles 

|:ReFo(^) = jyC^^3(C,M) < for C > 0, (5.28) 

with 

P2(C; X, -a) = - 2C(1 + 2x^a^) sin5 + (1 + 2x'^aY 
PsiC; X, 5) = S{3x + 0{6'))e + (2(2 - x) + 0{S'))e 

+ S{{10x^ + 4x- 24) + 0{S^))C + {8{x - 1){2 - x) + 0{6^)). 

(5.29) 

Indeed, P2{u;x,±a) > 0, since its discriminant, as a quadric in u, is < 0. 
Moreover P3(u;x,6) > 0, for < 5 < e{x) with e{x) sufficiently small, since 
in that case the coefficients of ...,u^ are positive (since x e (3/2,2)) and 
C >0. Thus the chosen path for w is a path of steepest descent. 
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w \ 




Figure 10: Second deformation of the previous last contours in Figure [9] for 
a = (T+. At the critical point Wc, the path of w has angles ±7r/4, ±37r/4, 
while the path of w leaves it with any angle in (7r/2, 37r/4). The value of q 
will be chosen during the analysis. For cr = (T_, one picks the mirror image 
of the figure above about the vertical line through Wc, with w and w also 
flipped; i.e., the dashed and solid lines are interchanged. 



Slope of the function —Fo(w) from Wc = a+ to q and Wc to e^'^^^oo. Consider 
the curves parameterized by 



— (JX 

w = Wc+ (e ± i)C, for e = ±1 and C > 0. (5.30) 

2 — X 



One verifies 



where, using the curve relation f l5.17p . 

PiiC, x) = ee + 2e + eC + -(x - 1)(2 - x) (5.32) 

X X 

and 



l^-«l' = 2(^) P|(C;^,-^)>o, 

|^-/3|' = 2(^]'p|(C;x,a)>0 



(5.33) 



with 

PKC; X, -a) := + eCil - 2x^a^) - 2xV^ + x - 1, (5.34) 
showing at once the denominator of (15.311) is > 0. 
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For £ = 1 , the polynomial P| in the numerator of (15.311) is > for ( > 0, 
because its coefficients are all > in the range 2 > x > 3/2 and for ^ > 0. 
Therefore the derivative (15.311) is < for C > 0. 

For e = — 1 , the polynomial P^{(;x) will be < for large enough ^ > 0. 
However in the range 2 > x > 3/2, 



- X + 6x - 4 > and < 



4(x-l)(2-x) 4 



—x"^ + 6a; — 4 



< 



11^ 



and thus for < C < Co with 



Co 



4(x- l)(2-a;) 
— + 6x — 4 



(5.35) 



(5.36) 



the cubic above is strictly positive 



e=-l X V — + 6x — 4 



> 0. 

(5.37) 

This is the reason why for e = — 1 we bend the path at q to be horizontal, 
with q set to be equal to 



— CTX 

q := Wc + Co- {±i - 1). 

2 — X 



(5.38) 



Slope of the function —Fo{w) from q to q — oo. Consider the horizontal line 
given by 

— CTX 

w = q-C-^ , C>0. (5.39) 

2 — X 

Then, in the range —a > and x G (3/2, 2), 



d ^ , , ax 



-ax 



(2-x)2(x2-6x + 4) 



xa^P6{(;x) 
\w — a\'^\w — 



< 
(5.40) 



with 



P6(C; x) = {2- xfQ2{C, x) + (2 - x)C'P7(a:) + C'Q2(C; x) 



(5.41) 



Indeed, (52 (Ci x) and Q2{Ci x) are quadratic polynomials in C, with coefficients 
polynomial in x and P7(x) is a seventh degree polynomial in x. All three 
coefficients of Q2{C',x) are > for 3/2 < x < 2, while P^lx) > also as 
long as 2 > X > 3/2. The coefficients of C° and C^ of Q2(C;^) are > for 
2 > X > 1.70, which moreover has a positive minimum in C for 2 > x > 1.70, 
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thus proving the assertion for 2 > x > 1.7. A httle numerics in fact shows 
we can remove the restriction 2 > x > 1.7 and deduce the inequahty for 
2 > X > 3/2. 

Thus we have also shown that the chosen path for w is of steep descent. 
Thus the steep descent method can be apphed along these curves where the 
maximum of ReFo(w), — ReFo(w) occur at the saddle point Wc- The main 
contribution comes from the integration over a ^-neighborhood of the critical 
point Wc for both w and w. For small 6, the error made is of order e"'^™' with 
/i ~ 5^. Let us therefore choose S = K'^m'^^^m'^ with any 7 E (0, 1/20) fixed 
(i.e., m~^/^ 3> 5 S> m~^/^). Then, the only non-vanishing contribution in the 
m — > cxD limit is given by the integrations with \w — wd < (5, \w — wd < 5. 
In these small neighborhoods, we can apply series expansions f l5.26p . After 
the change of variables 



z := Km 
we finally get for a = a±, 

1 

73, 



1/4, 



w 



Wr. 



(5.42) 



{2m) 



dz 



dz 



(27ri) 



dz 



2, — z e-2''/4+eiz^/2--«i5+i?i 
1 „z4/4+6l222/2-y22+i?2 



(5.43) 



dz — 

y z 



where Q{i) = exp (F2(wc, O-^m}^"^ + F^iiwc, Vi)m}^^ + 0{m~^^^)^ is the conju- 
gation given in (15. 5p . and where the Ri are error terms, to be discussed later. 
Note the involution 6*1 ^ — ^ "^^2, z —z between the two integrals 
on the right hand side of (15.431) . which also respect the integration paths. 
The error terms Ri include the following local contributions: 

(a) 0{'m~^/^) of (15. 7p (uniform for Vi in a bounded set), 

(b) 0{5) = (9(m~^/^) from ^4(1/;) in (15.260 (uniform for 9i in a bounded set), 

(c) 0{m5^) = 0{m^^"'^^^'^^^), which is the corrections in the series expansions 
of Fo{w) of order higher than 4, see (I5.26p . 

Indeed, (b) is immediate since F4 is linear (see (15. 8p ). To see that (c) holds, 
we need to control the fifth derivative of Fq at Wc- We have 



max 

\w — Wc\<S 



4 max 

\w — Wc\<S 



{w — a)^ {w — (3)^ 



< 



16 



(5.44) 



for m large enough. 
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Finally, taking the m — >■ cxd limit to (15.431) the error terms vanishes and at 
the same time the integrals extend to infinity. The fact that z is not exactly 

is irrelevant, since the result is identical as soon as the direction has an an- 
gle strictly smaller than 7r/4 to the imaginary axis. Similarly one can deform 
the 5-path as depicted in Figure 4. This ends the proof of Proposition 15.11 
and thus also of Theorem 11.51 □ 

Remark 5.2. For future use, we point out that the elliptic curve E has three 
points above x = oo, only one of which is real, namely 

2-1/6 

a = — -TTT- + . . . for X — > oo. (5.45) 

At this point at infinity, one has, using the estimate (15.451) . (3 — a = 
hm^^oo = -2^/^ and assuming T = ^ - ^ = 0, also a + f3 = 0. 
This implies that P = —a = —2^^^. Note how this contrasts with (a, x) = 
(0,3/2), {X,T) = (0,0), in which case a = (3 = 0. To summarize, near the 
real points on S, namely near x = 3/2 and x = oo, one has the following 
leading terms (set 7 := ^^)- 

(^,x)~(7(x-^)^0, («,/3)~0, (X,T)~0, !f^^(x-^)^ 

(^,^)~ (^^,ooj, (a,/3)~(2i-2^), (X,T)~(-2i,0), 

(5.46) 



6 Limit to the quintic kernel 

In this section we present a conjecture concerning the process that will occur 
in the situation illustrated in Figure [5l In Theorem 11.51 and, in particular, 
in formula (I5.25p . it was observed that when a,f3 ^ (and only then), the 
tips of the cusps (r, ~ (iTm^/^Xm^/^) tend to the same point and that 
Wc — Wi = ^'^2-3^^'' ~^ cube root of Fq{w) turns into a quartic root. 

This also means that the starting and end points a and b for the wanderers 
tend to coincide and that the line connecting both points becomes vertical 
and tangent to the ellipse, as described in Figure [H This corresponds to 
the first situation in (15.461) . We now pick the second situation in (I5.46p . for 
which the cube root of Fq{w) also turns into a quartic root. However, this 
forces the points a and 6 to be a bit beyond \/2n; this means in particular 
that a > b, which actually violates the condition a < 6 in Theorem 11.11 One 
can think of the passage from x = 3/2tox = ooasa transition process. 
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The most natural strategy would be to set a = 6 = v^2ri + ^2m and take 
m,n —>■ oo together. Then, under an appropriate scaling limit, we expect 
to get a process with a quintic kernel. Of course, there will be a parameter 
tuning regulating how close the two Airy fields come together. For example, 
if m = n, then we have to choose a = 2\/2n + 0{n~^^^) since the fiuctuations 
of the first n Brownian bridges alone live on the n~^^^ scale. 

Evidence in favor of Conjecture 11.71 To give some evidence to this 
conjecture, we present two pieces of rigorous mathematics, concerning the 
(one-time) kernel, with d < b, with time r absorbed into d, 



3/3,£ ('^\'^ / ui-b 

At'&;&) = -^/ 'i^ 'ZJl . . (6.1; 



(2«) ir.> Jr,i e-^'/s+SiS ^ - u; 



Proposition 6.1. The kernel K!^'', as in ^6. can he continued analytically 
to a new kernel K'^ , as in 116.^) . with same integrand as kernel 116.1]) . by 
moving d and b in the complex plane from their original position d < b to a 
new position b < d on the real line: 



{2TTi) 



^—u)-' ri+t2U) I I — — f 

e / ' \uj-aJ yio~b 



(ei; 6) = 7;r-;i / du I du ^^^^j^ ,. • (6-2 



UJ — LO 



integrated over contours VL and VL as in Figure [771 

Proof, b corresponds to the black dot and d to the white dot in Figures [H] 
and [121 the dashed line refers to the cu-integration and the solid line to the 
u)-integration. 

We noticed in Remark [52] that the elliptic curve introduced in (15.171) . 
contains another real point, namely one covering x = 00 for which {a, (3) = 
(2^/^, —2^^/^). This clearly violates the inequality d = am^^^ < b = /?m^/^, 
crucial for the derivation of the kernel (16.11) . 

• Keeping u fixed but arbitrary on the solid line, one sees that the dashed 
line of Figure [12] (a) can be deformed into the dashed lines of Figure [12] ((b) 
+ (c)). Then one notices that the (c)-contribution vanishes. Indeed, (i) if u) 
belongs to the solid line, outside the dashed circle, the w-integral vanishes, 
the integrand being holomorphic; (ii) if u belongs to the solid line, inside 
the dashed circle, one picks up a residue and thus the cu-integral equals 
_i_g(6-?i)'^j further integrated with regard to to, one obtains 



1 

2TTi 



du e^^^-^i)'^ = 0, 

olid circle of (c) 
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Figure 11: New contour and ^2, with the black dot = b and the 
white dot = a. The sohd hne VL refers to the integration of the a;-variable, 
while the dashed hne refers to the a;-integration. 




{d) / (e) " 

Figure 12: Representation of the deformation of the integration variables for 
the case b < a. All the contours are clockwise oriented, the black dot is 6, 
the white dot is a. The solid line refers to the integration of the cD-variable, 
while the dashed line for a;-integration. The contributions of (c) and (e) are 
exactly zero. 
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and thus the only contribution comes from (b). 

• At the next stage, picking an arbitrary uj G dashed contour (b), one 
deforms the solid contour (b) into the solid contours (d) + (e). In the same 
way, ii uo ^ dashed circle, the cD-integration contributes nothing, the inte- 
grand being holomorphic; if G dashed circle, the (D-integration contributes 
_i_g(6-€i)<^- further integrated with regard to a;, one obtains 

J dashed circle of (c) 

and thus the integration over the (d)-contour is the only contribution. Fi- 
nally, the solid and dashed contours of (d) can further be deformed into 
contours (f), thus leading to the contours of Figure [2], as the black dot b 
migrates to the right of the white dot d through the C-plane; this ends the 
proof of Proposition 16.11 □ 

Proposition 6.2. Consider the kernel K^{S,i]i2), as in ^6.2^) with d > b. 
Then, defining the scaling 

d = (2m)^/3 (1 + |em-2/5 + 

b = (2m)i/3 (-1 - i0m-2/5 + lr]m-^/^) (6.3) 
6 = -i2mf' (1 - l^m-^/^ - |(^. - TsO') rn'"') , 
one obtains, in the m — oo limit, the quintic kernel K^{^i,C,2), 



2-1/3^-2/15 f f e-^V3H2^ I [ uj-b 

lim 7^ I doj I duj 



m 



UJ 



{2TTif Jn Jq e-^'/3+«i- - \^cj - ay \u-b, 

where C and C are the paths defined in Figure The limit is uniform for 
9,ri,vi,V2 in a bounded set. 

Proof. We shall give the proof in the case of 77 = 0; the case 77 7^ is easy to 
implement. As in the case of the Pearcey process (see Theorem 1 1.5 1) , consider 
the scaling = Xm?^^, d = am^^^, b = (3m} /'^ and the change of integration 
variables uj = wm^^^, uj = wm}^^. Then, the kernel 06.21) becomes 

rn}/^ r r pmF{w)-niF(w) 

= 7—72/^^/^^ ^ (6-5) 

(27r2) Jn w-w 

with 

F{w) := -w^S + Xw + \n{w - f3) - \n{w - a). (6.6) 
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where Q and Q are the contours of Figure [TH with the black dot being j3 and 
the white dot a. Here, one imposes the property that F'{w) experiences a 
4-fold zero at some point Wc, with a, P, X real and a ^ Wc, P 7^ Wc, i.e. one 
requires all Ai = 0: 

-{w - a){w - P)F'{w) -{w- Wc)'^ =: Aqw^ + Aiw'^ + A2W + A3 

= (4we -a- (3)w^ + {a(3 - Qwl - X)w'^ (6.7) 
+ (4w^ + X{a + (3))w - Xa(3 + a - (3 - w^. 

The coefficients Aq = Ai = Q imply Wc = + P) and X = a/3 — 
and consequently A2 = — ^^(q; + /5)(5q;^ — 6aP + 5/?^) = 0, whose only real 
solution is given by a = — /? and thus Wc = and X = a/3 = — a^. For these 
values, one has A3 = —a{a^ — 2) = 0, implying a = — /5 = 2^/^, X = —2^/^ 
and Wc = 0. To summarize 

p = -l^l^ <wc = ^<a = and X = -2^^^. 

Note this solution corresponds precisely to the real point on the elliptic curve 
£, covering x = 00, as obtained on the second line of fl5.46p (see Remark [5^ . 
Since we have a quintic leading term ~ mw^, we make the change of variables 
w = m~^l^zci and w = m~^^^za. The precise coefficients are chosen in order 
to simplify the final formula. Indeed, with fl6.3p we obtain 

mF{w) = mF(0) + V2Z - Oz^/S + 1z^ + 0{z^mr'^l^, zm'^/^), (6.8) 

with the error uniform for 6,Vi in a bounded set. The prefactor in (16.51) . 
after the changes of variables, becomes m^^^m' 

-1/521/3(1 + C(m-2/5)), which 
cancels with the 2~^/^m"^/^^ in front of the l.h.s. of (16. 4p (as m — > 00). 
Except for the error terms, the result of the theorem would follow. 

What remains to be seen is that the higher order expansions in the series 
do not contribute. We do it by the steepest descent method as for the 
Pearcey case. Consider the curve parametrized hy w = e^^'^^^x. Then, for 
the function F, as in (16. 6p . , with a, P and X substituted, 

F{w) = --w^ - l?l^w + In (u; + 2^'^) - In (u; - 2^/=^) (6.9) 
3 

one checks 

d ^ ^, . (x^ cos(7r/5) + 1) P „ / X 

— ReF u; = -2 . \ , ^ [ \ < for all a; > 0, 6.10 

dx ^ ' x4 + 2x2cos(7r/5) + 1 ^ ' 

One then checks that along the dotted loop in Figure [HI ReF(w) — ReF(O) < 
for w 7^ 0. It is at once visible by superimposing the dashed contour of 



53 



Figure 13: Contourplot of the function Re(F(x + iy) — F{0)). The value is 
high in dark regions and low in light regions. 



Figure fm onto the contour plot, as in Figure fT3l Then along the curve given 

hy w = e^'^'^^^x, 

^{-ReF{w)) = -2 f (^^^;<^/5) + l) for all a: > 0, (6.11) 

dx^ ^ " a;4 + 2x2 cos (vr/S) + 1 ' ^ ' 

and along the solid loop in Figure [HI — ReF(w) + ReF(O) < for w 7^ 0. 
This shows that the curves have the steepest descent property. Thus, if we 
integrate (in w) around a 5-neighborhood of the origin, the error term will 
be only of order 0{e~^^) with /i ~ 5^. We choose 5 = m~^/^nrC for any 
7 G (0,2/35). Then, uniformly for 6,Vi,V2 in a bounded set, the error term 
0(2;'^m~2/5) = 0(m'''''~2/5) as m ^ 00, as ^ < 6m^^^/a < ni< /a. In the 
limit, the only part of the contour in Figure [TT| which contributes in the end 
are the 8 rays emanating from the origin, which one deforms so as to form 
consecutive angles 7r/5. This then yields the quintic kernel K^{9,ri;vi,V2) 
with the integration paths C and C of Figure El thus establishing Proposi- 
tion ESI □ 

Remark 6.3. It is an open problem to know whether the Fredholm deter- 
minant of the quintic kernel defines a probability. 
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